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Abstract 

We investigate charge transport in pentacene-graphene nanojunctions employing density 
functional theory (DFT) electronic structure calculations and the Landauer transport formal- 
ism. The results show that the unique electronic properties of graphene strongly influence the 
transport in the nanojunctions. In particular, edge states in graphene electrodes with zigzag 
termination result in additional transport channels close to the Fermi energy which deeply af- 
fects the conductance at small bias voltages. Investigating different linker groups as well as 
chemical substitution, we demonstrate how the transport properties are furthermore influenced 
by the molecule-lead coupling and the energy level lineup. 

Molecular junctions, consisting of single molecules that are chemically bound to electrodes, 
represent interesting systems to investigate and understand mechanisms of nonequilibrium trans- 
port processes at the nanoscale. Recent experimental-"— and theoretical—"— works have elu- 
cidated the conductance properties of molecular junctions and revealed a wealth of interesting 
transport phenomena such as Coulomb blockade, Kondo effects, negative differential resistance, 
vibronic effects and local heating, as well as switching and hysteresis. 

Most experiments on single-molecule conductance have employed metal electrodes. Carbon- 
based materials, such as carbon nanotubes or graphene, represent another promising class of ma- 
terials for electrodes. Studies of molecular contacts with carbon nanotubes have demonstrated a 
number of advantages—*^ as compared to metal electrodes, in particular, rigidity, mechanical sta- 
bility and thus a more precise control of the molecule-electrode contact geometries. Graphene, 
furthermore offers excellent electronic properties,— especially high electron mobility and, further- 
more, may facilitate optical addressability of the nanocontact for optoelectronic applications.— 
Electron transport through molecular films connecting graphene electrodes— and molecular junc- 
tions with few layer graphene electrodes— were recently reported. Furthermore, using electron 
beam lithography, all-carbon junctions consisting of a carbon chain between graphene electrodes 
were realized.— 

In this work, we study the conductance properties of molecule-graphene nanojunctions em- 
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ploying a combination of density functional theory (DFT) calculations and the Landauer transport 
formalism. Considering several different examples, we show that charge transport in molecule- 
graphene nanojunctions is strongly influenced by the edge structure of the graphene electrodes. In 
particular, for graphene electrodes with zigzag termination, edge states may result in additional 
transport channels close to the Fermi energy that affect the conductance at small voltages pro- 
foundly. These findings extend previous theoretical work on graphene-molecular junctions.— "— 




Figure 1: Molecule-graphene junctions with a covalently bound pentacene bridge investigated in 
this work, (a) Overall structure of a junction. Two types of molecular bridges without (b,d) and 
with (c,e) CN-groups are considered. Two different geometrical arrangements of the molecular 
bridge have been analyzed: twisted geometries (b,c), where the molecular bridge is weakly coupled 
to the leads through amide-binding groups and nearly flat geometries (d,e), where the bridge is 
strongly coupled via oxygen atoms to the leads. 

As prototypical systems, we consider molecular bridges based on pentacene as the central 
unit (see Fig. [T]). The extended ;r-electron system of pentacene facilitates electron transport and 
makes this system, in combination with its optical properties, a promising material for molecular 
devices. Pentacene-based systems have been used, for example, in thin-film organic electronic 
devices.—"^ In the graphene-molecule junctions depicted in FiglB? and c, the molecular bridge 
is bound covalently via amide linker groups to two graphene nanocontacts. This binding group 
has been used successfully for molecular nanojunctions with carbon nanotubes.— The coupling of 
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the molecule-graphene nanoj unction to leads and an external circuit are modeled by self energies 
and absorbing potentials (vide infra). To analyze the influence of graphene edge states on the 
conduction mechanism, we investigate both zigzag-terminated graphene nanocontacts, as depicted 
in FiglH and armchair termination. 

In addition to the basic pentacene molecular bridge, we consider CN- substituted molecular 
bridges (Fig{T}; and e). The electron- withdrawing CN-group allows one to shift the location of the 
molecular energy levels relative to the Fermi energy of graphene and thus to vary the conductance 
at small voltages. Due to the almost perpendicular geometry of the molecular bridge with respect 
to graphene, the molecular junctions depicted in FigsUJ? and c exhibit only weak coupling between 
the electronic ;r-orbitals of pentacene and graphene and thus low conduction. As an example for 
strongly coupled systems, we consider the systems depicted in Fig[l]l and e, where the binding 
to the nanocontacts is achieved via oxygen atoms forming oxepine rings resulting in an overall 
planar geometry and a strong electronic interaction between the central pentacene molecule and 
the graphene nanocontacts. 

The theoretical methodology used to study charge transport in the systems considered is based 
on a combination of DFT calculations to characterize the geometry and electronic structure of the 
systems and Landauer theory to simulate the transport properties. Closed shell DFT calculations 
using the B3LYP exchange correlation functional together with the 6-31G(d) basis set are used 
to model the different systems investigated in this work. To avoid dangling bonds, the graphene 
edges have been saturated using hydrogen atoms. During the optimization process, geometrical 
constraints have been imposed to prevent artificial distortions of the nanocontacts. All calculations 
were carried out using G09— and TURBOMOLE^ codes. Following previous work,— the Kohn- 
Sham matrix obtained for the relaxed junction is identified with the Hamiltonian of the system. To 
facilitate the transport calculations, the Hamiltonian matrix is partitioned into blocks that corre- 
spond to the molecular bridge left and right contacts /// and Hr and coupling between them 
Vim, Vmr and V„„, respectively, as described by Benesch et al.— . Direct coupling between the 
left and the right contacts is neglected. To model the effect that extended graphene contacts have in 
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the electron transport, we employ absorbing boundary conditions described by complex absorbing 
potentials —iWa (oc = /, r), which are added to the Hamiltonian in the graphene contacts. Thus, the 
overall Hamiltonian H reads 



H 



^Hi-iWi Vi,„ ^ 

^ml fiffi Vfiir 

y K,n H,-iWrj 



(1) 



The matrices describing the absorbing potentials are assumed to be diagonal, {Wa)ij = WaiSij. 
Following previous work,— a polynomial function was employed to model Wai- Accordingly, the 
elements of the matrices, which correspond to the rth atom of the nanocontact, are described by 
the formula 

Wai=a\xi-X^a\'^, (2) 

where x,- is the x-coordinate of the atom, assuming that the x axis is directed along the molecular 
bridge, and x^ is selected to be close to the point where the molecular bridge is connected to the 
nanocontact. To determine the empirical parameter a, test calculations over a wide range of values 
were carried out. On the basis of that, the value of a was chosen from a stable parameter range, 
where the results do not depend on a Kopf and Saalfrank— . 

The transport properties of the nanojunctions are described by the transmission function. On 
the basis of the model introduced above, the transmission function for an electron with energy E 
is given by the expression— 

t{E)=tr{^rG„^^,Gi}, (3) 
where the trace is taken over the molecular degrees of freedom. The transmission function close 
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to the Fermi energy Ef determines the conductance ^ = j^^t{Ef) at small bias voltages. 
In Eq. (??), G,n denotes the Green's function projected on the molecular bridge 

G„,{E) = [E-H,„-Li{E)-L,{E)]-\ (4) 
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The matrix r« = — 2Im{Ec<}, which characterizes the energy-dependent width of the molecular 
resonance states, is proportional to the imaginary part of the self energy La- The self energy 
Ea(£') describes the coupling of the molecular bridge to the graphene leads and is given by the 
expression 

i:aiE) = V,„a[E-Ha+iWa]-Wa,n. a = l,r. (5) 
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Figure 2: Transmission through a weakly coupled pentacene- graphene junction. The upper panel 
depicts the transmission function for the pentacene bridge (solid line). The eigenenergies of the 
bridge (dash-dotted line) are also shown. The lower panel demonstrates the influence of CN groups 
on the transmission. The differences between zigzag nanocontacts (dashed lines) and armchair 
nanocontacts (dotted line) are also shown. The inset highlights the edge induced peak with and 
without CN substitution. 



We first consider the amide-linked nanojunction with the basic pentacene bridge (see FigH}?). 
Fig. [2] (top) shows the transmission function of this junction in the range of energies close to 
the Fermi level. The transmission function exhibits a series of peak structures, which can be 
classified into two types: (i) isolated resonance peaks with a maximum transmission value of 
unity, in the following referred to as regular or molecular resonances, and (ii) a small double peak 
structure close to the Fermi energy (see also the inset in Fig©. The former are associated with 



6 



HOMO-3 HOMO-2 LUMO+4 




-3-2-10 1 2 3 

E-Ef[eV] 



Figure 3: Energy spectrum of the molecular junction depicted in Fig{T]:. The highlighted area 
contains 8 edge states (including HOMO-3 and HOMO-4) and one molecular state (LUMO-i-4), 
see text for details. Through CN-functionalization of the pentacene bridge, the molecular state 
(LUMO-i-4) is embedded into the family of edge states. The energies of the states depicted corre- 
spond to peaks at E = —0.2 eV and E = +0.2 eV in Fig 12] (bottom, dashed line). 

resonant electron transport through the discrete energy levels of the molecular bridge and thus 
appear at energies close to the eigenenergies of the molecular levels. The width of these peaks 
varies significantly, depending on the coupling of the specific molecular level to the graphene 
leads. The latter double peak structure close to the Fermi energy, on the other hand, cannot be 
related directly to electronic states of the molecular bridge. Furthermore, this peak structure does 
not correspond to direct tunneling between the graphene nanocontacts because this mechanism was 
explicitly excluded in our treatment (cf. Eq. ??). As discussed in detail below, we assign the double 
peak structure close to the Fermi energy to an additional transport channel caused by the edge states 
of the graphene nanocontacts (in the following referred to as edge induced transmission channel). 
It is worthwhile to mention that similar surface induced transport channels have also been found 
in carbon-nanotube-based molecular junctions-^^*^ as well as in molecular junctions with metal 
electrodes,—!^ where they are known as a metal-induced gap states. Due to the proximity of the 
peak structure to the Fermi energy, it dominates the conductance at small bias voltages. 

In the nanoj unction considered in the top panel of Figj2l the amplitude of the edge-induced 
transmission peak is rather small. The amplitude can be strongly enhanced if molecular levels 
exist close to the Fermi energy. This is demonstrated in the bottom panel of Figj2l which shows 
the transmission function for the CN-substituted pentacene bridge. The CN-substitution causes an 
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overall shift of the molecular levels to smaller energies. Asa consequence, the energy level asso- 
ciated with the LUMO of pentacene approaches the Fermi energy. As a result, the edge-induced 
transmission peaks acquire additional intensity that increases their maximal value to almost unity. 
Furthermore, the interaction between the molecular resonance state and the edge-induced trans- 
mission peaks results in a splitting of the double peak. 

To further analyze the edge-induced transmission channel, we have investigated the electronic 
states of the overall nanojunction, as represented by their molecular orbitals, in a range of energies 
close to the Fermi level (see Figj3]). These states can be classified in two groups: (a) states local- 
ized preferentially in the graphene contacts and (b) molecular states mostly localized on the bridge 
(see, e.g. LUMO-i-4 in Figl3]). A detailed analysis of the (a)-type states shows that many of them 
are localized at the edge of the graphene contacts (see e.g. HOMO-3, HOMO-2 in Figl3]). These 
states are closely related to the well-known edge states of zigzag-terminated graphene nanorib- 
bons.— Due to their localization at the molecule-graphene interface, they interact strongly with 
the molecular states and induce the additional transport channels. This can be seen, for exam- 
ple, in the spectrum depicted in Fig j3] bottom, where the near degenerate pair of states HOMO-2 
and HOMO-3 correspond to the induced transmission peak of the CN-substituted pentacene at 
E = —0.2 eV (see Fig. O. On the other hand, the molecular state LUMO-i-4 in Fig 13] represents a 
regular transmission channel and corresponds to the peak at E = +0.2 eV (see Fig. [2l). 

To unambiguously show that the edge-induced transmission channel appears as a consequence 
of the existence of edge states in zigzag-terminated graphene nanocontacts, we have compared 
these results with the ones obtained for molecular junctions based on armchair-terminated graphene 
nanocontacts. As is well-known, armchair edges of graphene do not possess edge states.— As a 
consequence, the transmission functions for both kinds of molecular junctions should differ, and 
the edge induced transmission channel should be missing in the armchair-terminated graphene 
nanocontact. As shown in Figl2] bottom, where the transmission function for a CN-substituted 
pentacene-graphene junction with armchair-terminated graphene nanocontacts is depicted, the 
edge- state-induced transmission channels are indeed not present for this type of molecular junc- 



8 



tion, whereas all of the regular peaks similar to those found in the zigzag-terminated graphene 
nanocontact appear. This result thus unambiguously demonstrates that these additional transmis- 
sion channels are induced by graphene edge states. 

In the molecular nanojunctions considered above, the overall coupling between the ;r-electron 
systems of pentacene and the graphene leads is rather small due to the twisted (almost perpendicu- 
lar) orientation of the molecular bridge with respect to the graphene nanocontacts (see FigHh and 
c). To study charge transport for stronger molecule-lead coupling, we have also investigated the 
molecular junctions d and f depicted in FigH] In these junctions, the pentacene bridge binds to 
graphene via oxygen atoms forming oxepine rings, which facilitates a nearly planar structure and 
therefore a stronger coupling. The transmission functions (see FigjH) again show a peak structure 
with regular peaks related to molecular resonances and edge -induced transmission channels close 
to the Fermi level. However, the peaks are significantly broader than those found in the amide- 
bound molecular junctions due to the stronger coupling between the graphene nanocontacts and 
the pentacene bridge. Furthermore, the overlap of some of the peaks results in values of the trans- 
mission function that are larger than unity. The larger coupling also induces a significantly higher 
transmission amplitude of the edge-induced transmission peaks close to the Fermi level. In the 
CN-substituted junction (see FigjH bottom), this is further enhanced by the presence of a molecu- 
lar resonance close to the Fermi level, a result already found for the amide-bound CN-substituted 
junction. Interference between this molecular resonance and the edge-induced structures results in 
a complex triple peak structure, a characteristic feature that can be observed in the transmission 
function at the Fermi energy (see FigJH bottom). 

In summary, we have investigated charge transport in molecule-graphene nanojunctions. The 
results obtained for different pentacene-based nanojunctions reveal the unique properties of graphene 
as material for electrodes. In particular, we have shown that charge transport is strongly influenced 
by the edge structure of the graphene electrodes. Graphene electrodes with zigzag termination have 
edge states that result in additional transport channels close to the Fermi level that strongly affect 
the conductance at small bias voltages. While the influence of graphene edge states on transport 
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Figure 4: Transmission through strongly coupled pentacene-graphene junctions. The upper and 
lower panel depict results without and with CN-substitution of the molecular bridge, respectively 
(see models d and e in Fig{T])- The arrows show how regular and induced channels interact in the 
strong coupling regime to form complex peak structures in the transmission function. 

in nanojunctions has been discussed before,—"^ the present results show, to the best of our knowl- 



edge for the first time, the importance of edge state effects in molecule-graphene nanojunctions in 



the resonant transport regime based on first principles calculations. Our results also show that the 



transport properties are strongly influenced by the effective coupling between the molecular bridge 



and the graphene nanocontacts as well as the position of the molecular energy levels, which can 



be controlled by the choice of different linker groups and appropriate functionalization of the pen- 



tacene molecular bridge, respectively. These findings may facilitate further experimental studies 



and possible future applications in nanoelectronic devices. 
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